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ABSTRACT 

The Galactic X-ray binary Cyg X-3 consists of a compact object wind-fed by a Wolf-Rayet 
(WR) star We determine its binary parameters and the mass-loss rate. Using the published 
infrared and X-ray data, we calculate the mass ratio and inclination. Then we fit the depen- 
dence of the mass-loss rate on the mass for WR stars of the WN type, and take into account 
a relation between the mass-loss rate and the binary period derivative, observed to be > in 
Cyg X-3. This allows us to solve for the masses, inclination and mass-loss rate in the system. 
Our obtained mass-loss rate is almost identical to that from two independent estimates and 
consistent with other ones, which strongly supports the validity of this solution. We obtain 
the WR and compact object masses as 10.3^2 gM0, 2.4;^j [M0, respectively. Thus, either a 
neutron star or a low-mass black hole can be present. However, the radio, infrared and X-ray 
properties of the system suggest that the compact object is a black hole. 
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1 INTRODUCTION 



In spite of its discovery already in 1966 jGiacconiet aD 11963), 
Cyg X-3 remains poorly understood. It is located a t a distance of 
d ^ 7-9 kpc in the Galactic p lane jPickevI 1 19831 : IPredehl et all 
120001 : iLing. Zhang & Tanell2009h . The donor appears to be a Wolf- 
Rayet (WR) star of the WN type Jvan Kerkwiik et al. 1992, 1996; 
Ivan KerkwiikI 1993l : lFender. Hanson & Poole vil999l) . However, the 
compact o bject may be either a neutron star (N S) or a black hole 
(BH, e.g.. iLommen et alJbOOSi: hereafter L05: Ivilhu etai]|2009l 
hereafter V09). 

Here, we self-consistently determine the binary parameters 
and the mass-loss rate, M, of the donor. We use published results for 
the radial velocities and a relationship between the period deriva- 
tive, P, and the total mass of the system. We close the equations by 
finding a tight correlation between M and the WR mass, Mwr. 



2 THE BINARY PARAMETERS AND MASS-LOSS RATE 

The best existing determination of the mass function of the compact 
object appears to be that of Hanson, Still & Fender (2000), who 
measured the radial velocity of the WR star as A^wr - 109+13 
km/s. Then, V09 studied orbital modulation of X-ray lines in the 
system using Chandra, in particular those originating close to the 
compact object. In their determination of the radial velocity of the 
compact object, Kc, they calculated the average based on 4, 8 and 
10 bins. We note that in the case of 8 bins, the two with the max- 
imum blueshift were removed. That case gave the radial velocity 
substantially less than the other ones, which reduced the resulting 



velocity average. Here, we use their values for the 4 and 10 bins 
with, in addition, the case of the Chandra observation of Obsid 
6601 from V09. The unweighted average of the three measure- 
ments gives then the compact-object velocity of Kc - 469 ± 88 
km/s. We have also repeated their fits using their data from fig. 7, 
obtaining values very similar to their, and consistent with the above 

average. Thus, we use this value of Kq below. 

With the binary period of P = 0.19969 d jSingh et alj[2002h . 
the mass functions are /wR = 2.13!i**Mo and/c = O.O27;j;°|^^M0. 
(The original value of Kc - 418 ± 123 km/s of V09 implies /wr = 
I-^I-o m'^O-) T^^^^ (for our value of /wr) implies the mass ratio 
and the inclination of 



0.23+ 



Mc _ A'wR 
A^WR Kc 

(/fwR + KcYKcP 



(0.69 ±0.12) 
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IOMq/ 
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(2) 



respectively (note some misprints in the corresponding equations of 
V09). Here, we have estimated the errors as corresponding to the 
extrema of the velocities. At Mwr = IOMq, / ^ 43°!". If ; > 30°, 
which is likely given the strong IR and X-ray orbital modulation, 
the above constraints imply Mwr ^ 42Mq and the compact-object 
mass of Mc < 9Mq. 

Based on their population synthesis and evolution calcula- 
tions, L05 claimed that the case of an NS can correspond only to 
the donor being a low-mass, < 1.5Mq, He star accreting at a high 
rate onto the NS through a Roche-lobe overflow. Such an He star 
has a negligible wind and the observed WR features originate then 
in an outflow centered on the NS. Then, both IR and X-ray spectral 
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features originate in the neutron-star surroundings. It is then not 
clear how to explain the orbital modulation of the continuum, the 
emission lines fluxes and the radial velocities. Specifically, both the 
lines originating in the WR-like outflow and around the NS should 
have simil ar dependencies on t he orbital phase, contrary to the ob- 
servations l lHanson et allllOOOl : V09). On the other hand, V09 self- 
consistently explain the radial velocity and flux modulations with 
the orbital phase for both the X-ray and near-lR features (lines and 
continuum) in the frame of a model consisting of a real WR star and 
a compact companion. In particular, if the minimum flux (X-ray, 
near-IR) occurs when the WR transits in front of the compact ob- 
ject, the ra dial velocity curves tracing the WR star (He: absorption; 
iHanson e t al. 2000) and the compact object (Fexxvi; V09) have the 
minima/maxima at the respective quadratures, i.e., are shifted by 
Q.25P and Q.15P, respectively, whereas the radial velocities of the 
IR Hen emission formed in the X-ray shadow cone behind the WR 
star, where the wind is photo-ionized by its EUV radiation, show 
a direct correlation with the light curve. Such a behavior appears 
impossible to explain in terms of the model in which all WR fea- 
tures arise in the matter re-ejected from the region centered at the 
compact object. We note that this model was proposed before the 
results of V09 were obtained. Ruling it out, the mass loss is via the 
standard WR stellar wind. 

We have then two additiona l constraints. The binary is de- 
tached jvan Kerkwiik et alj [T992h . and thus there is no Roche 
lobe overflow and accretion is from the wind. Most of the 
wind leaves the system carrying away its angular momentum. 
This leads to slowing do wn the binary rotation at a rate (e.g., 
iDavidsen & Ostrikedll974h . 



(3) 



P 2M 

P ^ ~m"' 

where M = Mwr + Mr. lOglev. Bell Bumell & Fended JioOlh criti- 
cally considered assumptions made in obtaining this relationship. 
Based on their discussion, there appear to be no major prob- 
lems with this method when applied to Cyg X-3. The most re- 
cent published estimate for P/P is (1.05 ± 0.04) x 10"'' yr"' 
jSin gh et al. 2002). It has been confirmed and improved by S. Ki- 
tamoto (private comm.) by including Suzaku observations, which 
yields PjP ^ (1.03 ± 0.02) x 10"'' yr"'. Then, M ^ IMPjP ^ 
19Mo(M/10-5Mo/yr). 

The second constr aint is of the mass- l oss ra te vs. mass re- 
lationship in WR stars. Ivan Kerkwiik et all ( |l996[) classified Cyg 
X-3 as the WN 4-8 stellar type, which is intermediate between 
the WNE type, almost H free, and WNL, with substantial fraction 
of H, on the basis of the I and K-band emission lines. Thus, we 
consider here WN stars regardless of their H content, though the 
M(Mwr) relationship we find is similar if we restrict our sample 
to the He fraction of K > 0.92. We use the sample of table 5 of 
Nugis & Lamers (2000, hereafter NLOO), containing 34 WN stars 
(31 of them with the type WN 4-8), and giving M corrected for 
wind clumping (unlike a number of previous studies), which we 
show in Fig. [T] The metallicity of all these stars is close to solar, 
Z = 0.0172-0.0176. In the sample, 29 stars have Mwr < 22M0, 
and 5 stars have Mwr = (4O-57)M0. We can see that the M(Mwr) 
relationship is rather steep for the former but it saturates for the 
latter. Since our considerations above point to a relatively moder- 
ate mass of the WR star, we fit only the former subsample. Us- 
ing a fitting method symmetric in Mwr and M (see appendix C of 
IZdziarski et^l2oT if), we obtain 
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Figure 1. The mass-loss rates corrected for wind clumping vs. the mass for 
WN stars in the sample of NLOO. The solid line shows the fit to all points 
with Mwr < 22Mo, and the dotted lines delineate the rms of these points. 
The red and blue crosses show the objects with reliable determinations of 
Mwr and M and those with only rough estimates (according to NLOO), 
respectively. We see no systematic difference between these two categories 
at MwR < 22Mo. 



(M) ^ 14.7Mo, 



2.93 ±0.38, M<22Mo, 



(4) 



Mn 



(1.9 ± 0.2) X 10" 



shown in Fig.[T] Here, the normalization is at the average of the fit- 
ted Mwr, (M). Equation Q for a fiducial M = lO"^M0/yr implies 
Mwr - (12 ± l)Mo. NLOO fitted the entire WN sample, includ- 
ing the massive stars, and thus obtained a much flatter dependence, 
see their equation (23). However, it is clear from Fig.[T]that this is 
a very inaccurate description of the actual dependence. Note that 
clumping reduces M for a given Mwr, and the above relationship 
is a factor of 2-3 below that of Langer (1989) or LOS. On the other 
hand, our fit predicts M higher by a factor of 3.6 than that given 
bv lVink & de Koteij ( |2005|) for the WN 8 star WR 40 (also present 
in the sample used here) with Mwr - 2OM0, used as the normal- 
ization of their theoretical dependence of M(Z) for WN stars. This, 
however, appears to be due to a mistake in their coping that value of 
M from Herald, Hillier & Schulte-Ladbeck (2001), whose value of 
3.1 X 10"^Mo/yr is almost the same as t hat of NLOO, and con sistent 
with our fit. We also note that although iHerald et ai] ( l200lh claim 
that their value of M for another WN 8 star, WR 16, is different 
from that of NLOO due to a different treatment of clumping, the dif- 
ference is entirely accounted for by the different values of the used 
distance. Summarizing, we have not found any valid objections to 
the results of NLOO for their WN sample. 

Equations l|3j and ^ can be solved for Mwr and M as func- 
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tions of q, equation ([T), 



M = ••- 



IPMo 
P{1 + q){M} 



2F 



(5) 
(6) 



Using the P/P of Kitamoto, these formulae give Mwr - 
10.3!]fjMo, Mc ^ 2.4+;;;5Mo, M ^ 6.5+|;^ X IQ-^Mo/yr, and 
i = 43°tg* [from equation (|2]l]. The lower limits on the masses 
and M and the upper limit on ; correspond to the lower limits on q, 
P and n and the upper limit on Mg, and vice versa. Using the orig- 
inal mass function of V09 changes these results only slightly. We 
note here that the fitting method finds the uncertainty of the aver- 
age dependence. Its spread, or rms, is much larger, and it is partly 
intrinsic and partly due to measurement errors. Without detailed 
knowledge of the errors, we can make a conservative assumption 
that the observed dispersion of M for a given Mwr is mostly in- 
trinsic, and due to the internal structure and/or Y affecting M. An 
estimate of the rms is AMo VA' - 1, where A' = 29 is the number of 
stars in the fitted sample. We then keep n fixed, which corresponds 
to moving up and down the best fit dependence within the observed 
spread, as illustrated by the dotted lines in Fig. [T] This way of 
estimating the M(Mwr) uncertainty yields somewhat wider con- 



fidence limits of Mwr ^ IO.S+^-^Mq, Mc 2.4+j |Mo, / = 43' 



f 17 
-12' 



M ^ 6.5^, , X 10 Mq/jt. Concluding, both ways of estimating the 
uncertainties allo w either a low-mass BH or an NS to be present. 

We note that ISingh et alj ( [20o3) and Kitamoto (priv. comm.) 
also fitted cubic ephemerides to their data, i.e., allowed for P ^ 0. 
They obtained rather small values of P/P ^ -60 ± 70 yr and 
-130 ± 80 yr, respectively. These ephemerides would thus predict 
that P changes its sign in ~ 10^ yr, which is obviously <c any evolu- 
tionary time scale of the wind. Thus, this model is inconsistent with 
the P being due to the loss of the angular momentum via wind, as 
assumed in equation l|3j- If such an ephemeris were shown to be 
correct, another explanation, e.g., a third body, would be required. 
However, a linear + sinusoidal ephemeris gives a very bad fit to 
those data. Moreover, a cubic ephemeris for the most recent data 
set yields no improvement of the fit over the parabolic one, and the 
same reduced^' (Kitamoto, priv. comm.). We thus do not consider 
this ephemeris any more. 

We then compare our results on M with other independent es- 
timates. The most detailed stu dy of the constraint from X- ray ab- 
sorption appears to be that of ISzostek & Zdziarskil ( 1200 8l) . They 
took into account wind clumping as well as the presence of a ho- 
mogeneous wind phase. They obtained results for the cases of a 
low and high mass of the WR star, 6.4Mq and 70Mq, respectively. 
Our range of Mwr is much closer to the lower one, and the depen- 
dence on Mwr is relatively modest. For the low-mass case, they 
obtained M (6-8) x 10"''MQ/yr, almost the same as our val- 
ues obtained here. A related (but independent) estimate is that of 
Zdziarski et al. (2012), who studied X-ray orbital modulation in 
Cyg X-3. Its orbital-phase dependence in hard X-rays, where the 
opacity is likely to be dominated by Thomson scattering, implies 
M ^ (7-8) X lO^'^Mp/yr for the terminal wind velocity o f 1500- 
1700 km/s (lFenderetal.lll999l ; lvanKerkwiiketai]|l996i) . Then, 
IWaltman et alj i 1 1 996h obtained M < 1 .Ox lO^^Mo/yr from the radio 
delays interpreted as being due to wind opacity. This is again fully 
consistent with our values. The above methods are independent of 

the distance to Cyg X-3. 

On the other hand, Ivan KerkwiikI ( Il993h and lOglev et"al] 

l l200lh obtained higher values of M by interpreting the observed 



IR fluxes as solely due to wind emission. Ivan KerkwiikI ( jT993h ob- 
tained M 4 x IO"^Mo/yr without taking into account clump- 
ing of the wind. Including it lowers M by a factor of several 
jSzostek & Zdziarskill2008l) and thus makes it consistent with our 
present results. Ogley et al. (2001) obtained M ^ (4-30) x 
lO^'MQ/yr using different assumptions about the wind, but also ne- 
glecting clumping. The lowest value is for a strongly ionized wind. 
Using that assumption and taking into account clumping can again 
reconcile their results with ours. Furthermore, the wind-emission 
calculations depend on the distance as M oc d''^, and the above au- 
thors assumed d - 10 kpc. \f d ^ 1 kpc, their values of M should 
be reduced by ~40 per cent. In addition, there can be a substantial 
IR contribution from the synchrotron emission of the jets, which 
would again reduce the calculated M. Summarizing, all the pub- 
lished estimates of M are compatible with our present results. 



3 DISCUSSION 

As we noted above, L05 claimed that only a BH can be present in 
the case of a WR donor, based on their population synthesis results. 
However, given uncertainties intrinsic to the population-synthesis 
method, this conclusion appears uncertain. Also, their specific ver- 
sion of this method suffers from a number of problems. One of 
them is their condition for accretion disc formation for an isotropic 
wind of Ergma & Yungelson ( 1998). However, if the wind veloc- 
ity at the Roche-lobe surface of the WR donor is less than the 
escape velocity from that surface, the wind will be strongly fo- 
cused towards the accretor, which facilitates accretion disc forma- 
tion. Numerical simulations for such a wind Roche-lobe overflow 
(Mohamed & Podsiadlowski 2012) indicate that this mass transfer 
mechanism may occur in systems where the wind acceleration zone 
lies close to or makes a significant fraction of the Roche-lobe ra- 
dius. This is very likely to be the case in Cyg X-3 given its short or- 
bital period. Thus, our conclusion above that an NS may be present 
in Cyg X-3 appears to be valid. 

On the other hand, there is a lot of other (though still cir- 
cumstantial) evidence for the compact object in Cyg X-3 being 
a BH. Its spectral states have spectra showing a general resem- 
blan ce to those of the canonical spectral states of BH binaries 
(iSzostek & Zdziarskil l2008l : IZdziarski. Misra & Gierhhski 2O10i). 
The existing difl^erences may be explained by scattering in the very 

The radio 



The existing difl^erences may be explained by scattering i 
strong stellar wind in Cyg X-3 jZdziarski et alJlToid) . 
flux in the hard spectral state is strongly c orrelated with that in soft 
X-rays, which correlation is very simi lar ( Gallo . Fender & PoolevI 
l2003l ; ISzostek. Zdziarski & McCollough,2008.) to that seen in BH 
X-ray binaries. The normalization of the radio flux in Cyg X-3 is 
somewhat above that for BH binaries, which excess can be due to 
strong absorption of soft X-rays in Cyg X-3. On the other hand, 
NS X-ray binaries in the analogous state have the r adio flux > 30 



times weaker on average jMigliari & Fenden 



20061). in stark con- 



trast to the case of Cyg X-3. As pointed out by Szostek & Zdziarskil 
(12008^ all NSs in confirmed high mass X-ray binaries are highly 
magnetized (which is related to them being young systems), and 
are usually X-ray pulsars. Even in cases with no detected X-ray 
pulsations, their bremsstrahlung type spectra bear no resemblance 
to Cyg X-3, which is also a young system. Then, V09 fitted the 
IR-flux orbital modulation and found it is compatible with a low 
inclination, i ~ 30°, but not with ; ~ 60°, implying the compact ob- 
jec t to be a BH. A l ow in clination, / ^ 30-40°, has also been found 
bv lZdziarski et alj ( 1201 2h as accounting for the orbital modulation 
of the X-ray continuum. 
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Following the above circumstantial evidence and our esti- 
mate of the compact object mass (Mc = 2.4^| [Mq) we ar- 
gue that Cyg X-3 hosts a low mass BH. The actual BH mass 
is most likely in the range (MNs,max-4.5MQ), where MNs.max is 
the still unknown maximum NS mass. The recent precise mass 
measurements of Gal actic pulsars require that MNs.max > 2Mq 
jPemorest et alj|201(l) . w hile dynamical BH mass estimat es im- 
pose that MNs,max < 5Mq jOzel et alfcOld l lFarr et aTll201lh . The- 
oretical calcul ations of den se nuclear matter indicate that Mns.hmx < 
2.4Mq (Lattimer & Prakash 2010 ). When we impose above limits 
on our estimate we obtain Mbh = 2.4^g _[MQ. This puts the BH in 
Cyg X-3 right into the mass gap, the d earth of compact objects in 
the (2-5)Mq mass range, first noted bv lBailvn e t al. ( 1998). 

The existence of the mass gap w as recently discussed in con- 
text of supernova explosion models teelczviiskietalj|2012h . For 
the rapidly developing explosions (within ~200 ms after the core 
bounce), the mass gap appears in simulations as observed for the 
Galactic BH X-ray binaries. However, there is a small addition 
of a few syst ems with compact obje cts in (2-3.5)Mq mass range 
(see fig. 1 of iBelczvnski et alj|2012l) . These have formed via ac- 
cretion onto an NS (originally formed in the supernova explo- 
sion) from a binary companion. These compact objects are most 
likely light BHs formed via accret ion induced collapse of the NS 
jvan den Heuvel & De Loordll973h . 

The mass gap was also recently interpreted as a potential ob- 
servational artefact caused by th e systematic uncertain ties in dy- 
namical BH mass measurements ( Kreidb erg et alj|2012t) . No mass 
gap scenario corresponds to delayed supernova models (develop- 
ing on a timescale ~500-10 00 ms) that produce co mpact objects 
within (2~5)Mq mass range i Belczviiski et al.|[2012l) . In this case, 
the light BH in Cyg X-3 most likely have formed in su pernova ex- 
plosio n from an Mzams = (20~40)Mq progenitor star jFrver et al] 
I2OI2I) . As long as there is no consensus on the mass gap existence, 
both formation scenarios, either via accretion induced collapse or 
in the delayed supernova explosion, are to be considered. 

The Eddington luminosity for He accretion is Le - 2.5 x 
10''*(Mc/Mq) erg s"'. The absorption corrected, phase-averaged, 
bolometric fluxes in th e hard and soft spectral s tate (assuming 
the low-mass m odel of Szostek & Zdziarskill2008l : the hard-state 
flux corrected in ISzostek et alj|2008l) are 8.5 x 10 erg cm s 
and 2.7 x 10"** erg cm"^ s"', respectively. These are orbital-phase 
averages; they should be increased by a factor 1.5 and 2 in 
the hard and soft stat e, respectively, to obtain the intrinsic values 
jZdziarski et al.ll2012l) . Assuming d - 7 kpc, we obtain the bolo- 
metric hard and soft state luminosities of ^ 7 x 10^' erg s"' and 
- 3 X 10^** erg s"', respectively. At M = 4Mq, they correspond 
to ^ 0.07Le and - 0.3Le, respectively. These Eddington ratios 
are wi thin the ranges observed for black -hole binaries in these two 
states toone, GierUiiski & Kubotall2007l) . 



4 CONCLUSIONS 

We have have determined the masses of the components, the incli- 
nation and the donor mass-lo ss rate in the Cyg X-3 binary. We have 
used the Kq measurement of iHanson et al.l ( I2OO0I) and revised the 
determination of A'wr of V09 based on their data. We then assumed 
that M is related to P and the total mass via the orbital angular 
momentum being lost by the wind. We have also substantially re- 
vised existing dependencies of M(Mwr) for WN stars by using the 
dumping-corrected results of NLOO and fitting them for moderate 
masses only. The resulting four equations are solved for Mwr, Mc, 



i and M. The obtained M is in a very good agreement with two 
independent calculations, and are consistent with all other deter- 
minations, provided those based on the IR fluxes are corrected for 
clumping. 

Our results for the (conservative) assumption that the ob- 
served spread of M(Mwr) is intrinsic are Mwr - 10.3;^^gMQ, 
Mc ^ 2.4^2 |Mo, = 43°*\l, M ^ 6.51^-J x lO^^Mo/yr. We ar- 
gue that although an NS is allowed to be present, a BH is more 
likely given the radio, IR and X-ray properties of Cyg X-3. 
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